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We analyze the constrained minimal supersymmetric standard model with /i > supplemented 
by a generalized 'asymptotic' Yukawa coupling quasi-unification condition, which allows an accept- 
able 6-quark mass. We impose constraints from the cold dark matter abundance in the universe, B 
physics, and the mass of the lightest neutral CP-even Higgs boson. We find that, in contrast to 
previous results with a more restrictive Yukawa quasi-unification condition, the lightest neutralino 
X can act as a cold dark matter candidate in a relatively wide parameter range. In this range, the 
lightest neutralino relic abundance is drastically reduced mainly by stau-antistau coannihilations 
and, thus, the upper bound on this abundance from cold dark matter considerations becomes com- 
patible with the recent data on the branching ratio of Bs — > fi'^ . Also, ~ (125 — 126) GeV, 
favored by LHC, can be easily accommodated. The mass of x, though, comes out large (~ 1 TeV). 

PACS numbers: 12.10.Kt, 12.60.Jv, 95.35.+d 



I. INTRODUCTION 

The recently announced experimental data on the mass 
of the standard model (SM)-like Higgs boson jH-lll as well 
as the branching ratio BR(i?s — > iJi~) of the process 
Bg — >■ H'^ fJ-~ 1^] in conjunction with cold dark matter 
(CDM) considerations Q put under considerable stress 
[6| the parameter space of the constrained minimal super- 
symmetric standard model (CMSSM) [7l-[l0j. Let us re- 
call that the CMSSM is a highly predictive version of the 
minimal supersymmetric standard model (MSSM) based 
on universal boundary conditions for the soft supersym- 
metry (SUSY) breaking parameters. The free parameters 
of the CMSSM are 



sign/x, tan /?, M- 



1/27 ""^0, 



and An 



(1) 



where sign/i is the sign of /i, the mass parameter mix- 
ing the electroweak Higgs superfields H2 and Hi of the 
MSSM which couple to the up- and down-type quarks 
respectively, tan /3 is the ratio of the vacuum expectation 
values (VEVs) of H2 and Hi , and the remaining symbols 
above denote the common gaugino mass, the common 
scalar mass, and the common trilinear scalar coupling 
constant, respectively, defined at the grand unified the- 
ory (GUT) scale Mqut determined by the unification of 
the gauge coupling constants. 

It would be interesting to investigate the consequences 
of these experimental findings for even more restricted 
versions of the CMSSM which can emerge by embedding 
it in a SUSY GUT model with a gau ge g roup containing 
SU{A)c and SU{2)ii. This can lead [11[ to 'asymptotic' 
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Yukawa unification (YU) |12| . i.e. the exact unification 
of the third generation Yukawa coupling constants (of 
the top [bottom] quark ht [hb] and the tau lepton hr) 
at Mqut- The conditions for this to hold are that the 
electroweak Higgs superfields H2, Hi as well as the third 
generation right handed quark superfields form 5/7(2)^ 
doublets, the third generation quark and lepton SU{2)l 
doublets [singlets] form a SU{4)c 4-plet [4-plet], and the 
Higgs doublet Hi which couples to them is a SU{A)c 
singlet. The simplest GUT gauge group which contains 
both SU{4)c and SU{2)r is the Pati-Salam (PS) group 
Gps = 5C/(4), X SU{2)l X SU{2)r 

It is well-known that, given the experimental values of 
the top-quark and tau-lepton masses (which, combined 
with YU, naturally restrict tan/3 to large values), the 
CMSSM supplemented by the assumption of YU yields 
unacceptable values of the 6-quark mass for both signs 
of /i. This is due to the generation of sizable SUSY correc- 
tions 15] to TTifc (about 20%), which arise from sbottom- 
gluino (mainly) and stop-chargino loops [l^ [l^ and have 
the same sign as /i - with the standard sign convention of 
Ref. • The predicted tree-level mi,{Mz), which turns 
out to be close to the upper edge of its 95% confidence 
level (c.l.) experimental range, receives, for fi > [fi < 0], 
large positive [negative] corrections which drive it well 
above [a little below] the allowed range. Consequently, 
for both signs of /z, YU leads to an unacceptable mb{Mz) 
with the fj, < case being much less disfavored. 

In Ref. [13 - see also Refs. [Hil -, concrete SUSY 
GUT models based on Gps are constructed which nat- 
urally yield a moderate deviation from exact YU and, 
thus, can allow acceptable values of the 6-quark mass 
for both signs of /i within the CMSSM. In par ticular, 
the Higgs sector of the simplest PS model [Tsl , [T^ is 
extended so that H2 and Hi are not exclusively con- 
tained in a 5[/(4)c singlet, SU{2)l x SU{2)r bidoublet 
superfield, but receive subdominant contributions from 
another bidoublet too which belongs to the adjoint rep- 
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resentation of SU{4:)c- As a consequence, a modest vio- 
lation of YU is naturally obtained, which can allow ac- 
ceptable values of the 6-quark mass even with universal 
boundary conditions. This approach is an alternative 
to the usual strategy f23l - l26j according to which YU is 
preserved, but the universal boundary conditions of the 
CMSSM are abandoned. 

We wiU focus, as usually [l8l [20l - l22l | . on the fi > case 
since /i < is strongly disfavored by the constraint aris- 
ing from the deviation Sa^ of the measured value of the 
muon anomalous magnetic moment from its predicted 
value a™ in the SM. Indeed, ^ < is defended [131 only 
at 3 — cr by the calculation of a™ based on the r-decay 
data, whereas there is a stronger and stronger tendency 
[28l [29} at present to prefer the e''"e~-annihilation data 
for the calculation of a^^, which favor the /z > regime. 
Note that the results of Ref. [13], where it is claimed 
that the mismatch between the r- and e~'"e~-based cal- 
culations is alleviated, disfavor ^ < even more strongly. 

The representation used for the Higgs superfield which 
mixes the SU{2)l doublets contained in the SU{4:)c sin- 
glet and non-singlet Higgs bidoublets plays a crucial 
role in the proposal of Ref. [13]. As argued there, this 
Higgs superfield can be either a triplet or a singlet under 
SU{2)ii. In particular, it was shown that extending the 
PS model so as to include a pair of S'C/(2)fl-triplet and/or 
a pair of S't/(2)i^-singlet Higgs superfields belonging to 
the adjoint representation of SU (4)c can lead to a sizable 
violation of YU. However, in the past, we mainly focused 
[isl [20l - [2^ on the minimal extension of the PS model re- 
sulting from the inclusion of just a pair of SU (2) n-triplet 
superfields since this was enough to generate an adequate 
violation of YU ensuring, at the same time, a SUSY spec- 
trum which leads to successful radiative electroweak sym- 
metry breaking and a neutralino lightest SUSY particle 
(LSP) in a large fraction of the parametric space. The re- 
sulting asymptotic Yukawa quasi-unification conditions, 
which replaced the exact YU conditions, depend only on 
one new complex parameter (c) which was considered 
for simplicity real. It is also remarkable that this model 
predicts [2l|, ll^j values for the mass mhof the CP-even 
Higgs boson h close to those discovered [ll-Q by the Large 
Hadron Collider (LHC) and supports new successful ver- 
sions [3lj of the F-term hybrid inflation based solely on 
renormalizable superpotential terms. 

However, it has been recently recognized 22] that the 
lightest neutralino x cannot act as a CDM candidate in 
this model. This is because the upper bound on the light- 
est neutralino relic density from CDM considerations, 
although this density is strongly reduced by neutralino- 
stau coannihilations, yields a very stringent upper bound 
on the mass of the lightest neutralino m^, which is incom- 
patible with the lower bound on m-^ from the data [32] on 
BR(_Bs — > This result is further strengthened by 

the recent measurements [3] on BR(i3s — >■ which 
reduce the previous upper bound on this branching ratio 
and, thus, enhance even further the resulting lower bound 
on m^. The main reason for this negative result is that 



tan /3 remains large and, thus, the SUSY contribution to 
BK{Bs — >■ which originates [13, H^] from neutral 

Higgs bosons in chargino-, H"^-, and W^^-mediated pen- 
guins and behaves as tan^ f3/m\, turns out to be too 
large (m^ is the mass of CP-odd Higgs boson) . 

Therefore, it would be interesting to check if, in the 
framework of the CMSSM and consistently with the GUT 
models of Ref. ]1§] , we can revitalize the candidacy of x 
as a CDM particle, circumventing the constraint from 
BR(_Bs — and, at the same time, obtaining ex- 
perimentally acceptable rrihS. A key-role in our present 
investigation is the inclusion of both pairs of SU{2)ii- 
triplet and singlet Higgs superfields. This allows for a 
more general version of the Yukawa quasi-unification con- 
ditions - already extracted in Ref. [l^] - which now de- 
pend on one real and two complex parameters. This lib- 
erates the third generation Yukawa coupling constants 
from the stringent constraint hb/ht + hr/ht = 2 ob- 
tained in the monoparametric case and, thus, can ac- 
commodate more general values of the ratios hi/hj with 
I, j = t,b,T, which are expected, of course, to be of order 
unity for natural values of the model parameters. This 
allows for lower tan /3's and, consequently, the extracted 
BR(i3s — >• can be reduced to an acceptable level 

compatible with the CDM requirement. The allowed pa- 
rameter space of the model is then mainly determined by 
the interplay of the constraints from BR(i?s — > /i^/i") 
and CDM and the recently announced results of LHC on 
the Higgs mass m/j. 

We first review the salient features of the PS GUT 
model in Sec. |lT] and exhibit the cosmological and phe- 
nomenological requirements that we consider in our in- 
vestigation in Sec. IIIII We then find the resulting restric- 
tions on the parameter space of our model and test the 
perspective of direct neutralino detectability in Sec. IIVI 
Finally, we summarize our conclusions in Sec. [V] 



II. VIOLATING YU WITHIN A SUSY PS 
MODEL 

The starting point of our construction is the SUSY 
GUT model presented in Ref. [13]. It is based on Gps, 
which, as already mentioned, is the simplest gauge group 
that can lead to YU. The representations and transforma- 
tions under Gps of the various matter and Higgs super- 
fields of the model are presented in Table|T]([/c G SU{A)c, 
Ul e SU{2)l, Ur e SU{2)r and t ,t, and * stand for 
the transpose, the hermitian conjugate, and the complex 
conjugate of a matrix respectively). The model also pos- 
sesses a Peccei-Quinn (PQ) symmetry, a U{1) R symme- 
try, and a discrete matter parity symmetry with the 
charges of the superfields under these extra global sym- 
metries also shown in Table HI The matter superfields are 
Fi and [i — 1,2,3), while Hi and H2 belong to the 
superfield li. So, as one can easily see, all the require- 
ments [ll| for exact YU are fulfilled. The breaking of 
Gps down to the SM gauge group Gsm is achieved by 
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TABLE I: Superfield Content of the Model 



Super- 


Represen- 


Transfor- 




Global 




fields 


tations 


mations 


Symmetries 




under Gp; 


3 under Gps 
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PQ 




Matter Fields 




(4,2,1) 


FUluJ 


1/2 


-f 


1 


Fr 


(4,1,2) 


U*UiFr 


1/2 





-1 


Ifiggs Fields 


H" 


(4,1,2) 













H" 


(4,1,2) 


H^UluJ 











S 


(1,1,1) 


S 


1 








G 


(6,1,1) 


UcGUj 


1 








h 


(1,2,2) 







1 





N 


(1,1,1) 


N 


1/2 


-f 





N 


(1,1,1) 


N 





1 





Extra Higgs Fields 


li' 


(15,2,2) 


U*ULlh'u1u] 





f 





h' 


(15,2,2) 


UMlIi'uIu^ 


i 


-i 





^ 


(15,1,3) 















(15,1,3) 




f 








4>' 


(15,1,1) 













0' 


(15,1,1) 




i 









the superheavy VEVs (~ A/qut) of the right handed 
neutrino type components {vj^, i>jj) of a conjugate pair 
of Higgs superfields H'^, H'^. The model also contains 
a gauge singlet 5" which triggers the breaking of Gps, 
a SU(4)^ 6-plet G which gives [13] masses to the right 
handed down quark type components of H^, H'^, and a 
pair of gauge singlets N, N for solving [s^ the n problem 
of the MSSM via a PQ symmetry. 

In order to allow for a sizable violation of YU, we ex- 
tend the model by including three extra pairs of Higgs 
superfields lh',]h , (j),(f>, and where the barred su- 

perfields are included in order to give superheavy masses 
to the unbarred superfields. These extra Higgs super- 
fields together with their transformation properties and 
charges are also included in Table HI The superfield ]h' 
belongs to the (15,2,2) representation of SU{4:)c which 
is the only representation, besides (1,2,2), that can cou- 
ple to the fermions. On the other hand, (j> and (j>' acquire 
superheavy VEVs of order A/qut after the breaking of 
Gps to GsM- Their couplings with Jh and Jh naturally 
generate a SU(2)^- and SU (4) ^-violating mixing of the 
SU{2)l doublets in Ih and Ih' leading, thereby, to a siz- 
able violation of YU. 

More explicitly, the part of the superpotential which is 
relevant for the breaking of Gps to Gsm is given by 

-S{l3(b''+l3'(b'^) + {\^ + X'^')H^H', (2) 

where the mass parameters Af, to, and m' are of order 
-^GUT, and K, /3, /?', A, and A' are dimensionless param- 
eters with Af , TO, to', k, a, A' > by field redefinitions. 



For simplicity, we take /? > and /?' > (the parameters 
are normalized so that they correspond to the couplings 
between the SM singlet components of the superfields) . 
The scalar potential obtained from Wh is given by 



'J.'2 I 



iKiH^H" - Nr) - I3(j3^ - 13' (j)' 
+ |K5 + A^ + AV|'(|i7"P + |i5-^|^ 



\2liS<t>-m(t>\ 
\m(l) + XH'=H'= 
D — terms. 



\2l3'S(f>' - m'(f>'\ 
f \m'(j)' + X'H'^H'^ 



(3) 



where the complex scalar fields which belong to the SM 
singlet components of the superfields are denoted by the 
same symbols as the corresponding superfields. Vanish- 
ing of the D-terms yields H'= * = e'^H" {H", H" lie in the 
v"^, Djj direction). We restrict ourselves to the direction 
with -d — which contains the SUSY vacua (see below). 
Performing appropriate R and gauge transformations, we 
bring H'^, H'^ and S to the positive real axis. 

/^From the potential in Eq. ([3]) , we find that the SUSY 
vacuum lies at 



and 



where 



{S) = (^) = (^') = 0, 



f-V 

W0 = - 



= ^(l-^/TTl^ 

TO ^ to' 



with 



M2 

K 



I XI2 



< 1/4. 



(4a) 
(4b) 

(4c) 

(4d) 

(4e) 
(4f) 

(4g) 



The structure of (0) and (cf)') with respect to (w.r.t.) Gps 
is shown in Eqs. (l4b| and (l4c)) . where 



T, 



15 



1 



2^/3 



diag(l,l,l,-3). 



0-3 = diag (1,-1), and ctq = diag (1, 1) . 



(5a) 
(5b) 



The part of the superpotential which is responsible for 
the mixing of the doublets in Ih and Ih' is 



W,„ = Mhhlh' + X3(l)]h']h + Xx(j}'h!h, 



(6) 



where the mass parameter My^ is of order A^qut (made 
real and positive by field rephasing) and A3, Ai are 
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dimensionless complex coupling constants. Note that 
the two last terms in the right hand side (RHS) of 
Eq. ([6]) overshadow the corresponding ones from the 
non-renormalizable S'J7(2)fl-triplet and singlet couplings 
originating from the symbolic coupling H'^H'^Ii' Ih (see 
Ref. [3). Defining properly [l^, [2^ the relevant cou- 
plings in the RHS of Eq. ([6]) , we obtain the mass terms 



(7) 



where e is the 2x2 antisymmetric matrix with £12 = 1, 
the ellipsis includes color non-singlet components of the 
superfields, and the complex dimensionless parameters 
ai and a2 are given by 



1 



ai 



0-2 



1 



-f Ai-y^) , (8a) 
{\zv^ + Aiw;) • (8b) 



It is obvious from Eq. ^ that we obtain two pairs of 
superheavy doublets with mass My^: 



]h'i, H2 and H[, lh'2, 



where 



Ih' + arlhr 



1,2 



(9a) 



(9b) 



\ar\ 



(no summation over the repeated index r is implied) . The 
electroweak doublets Hr^ which remain massless at the 
GUT scale, are orthogonal to the H'^. directions: 



-allhL + Ihr 



(10) 



Solving Eqs. ^h\\ and ([TU| w.r.t. Ih^ and Ih'^, we obtain 

Hr + alH' I —arHr + Hi 

^=^= and Jh^ = 



Ihr = 



(11) 



The superheavy doublets H'r must have zero VEVs, 
which gives 



(.Hr 



and (li^) = 



-ar(Hr 



ar 



(12) 



The Yukawa interactions of the third family of fermions 
are described by the superpotential terms 



(13) 



where the factor of two is incorporated in the second term 
in the RHS of this equation in order to make j/33 directly 
comparable to 2/33, since the doublets in ]h' are propor- 
tional to T^^, which is normalized so that the trace of 



its square equals unity. From Eqs. p2)) and ([13]) and us- 
ing the fact that Ih' is proportional to T^^ in the SU{4:)c 
space, we can readily derive the masses of the third gen- 
eration fermions: 



mt 



nib 



1 - pa2/\/3 



(1 + KP) 

1 - pai/y/3 



2/33^2 



(l+|«l|2)^ 

1 -I- V^pai 



(l + l«iP) 



(14a) 
(14b) 
(14c) 



where p = 2/33/2/33 can be made real and positive by read- 
justing the phases of Ih, Ih' and Vr = (Hr). The third 
generation Yukawa coupling constants {ht,hb, and hr) 
must then obey the following set of generalized asymp- 
totic Yukawa quasi-unification conditions: 



htiMouT) : hb{MGVT) ■ hr{MG\JT) 



1 - P0L2l\fi 



1 - pai/VS 



Vi + KP 



1 + \fipa\ 



(15) 



These conditions depend on two complex (ai, 02) and 
one real and positive (p) parameter. For natural values 
of p, ai, and 012, i-C- for values of these parameters which 
are of order unity and do not lead to unnaturally small 
numerators in the RHS of Eq. (fTS]) . we expect all the 
ratios hi/hj with z,j = i, 6, r to be of order unity. So, 
exact YU is naturally broken, but not completely lost 
since the ratios of the Yukawa coupling constants remain 
of order unity restricting, thereby, tan/3 to rather large 
values. On the other hand, these ratios do not have to 
obey any exact relation among themselves as in the previ- 
ously studied p^l - l22j monoparametric case. This gives us 
an extra freedom which allows us to satisfy all the phe- 
nomenological and cosmological requirements with the 
lightest neutralino contributing to CDM. 



III. COSMOLOGICAL AND 
PHENOMENOLOGICAL CONSTRAINTS 

The two-loop renormalization group equations for the 
Yukawa and the gauge coupling constants and the one- 
loop ones for the soft SUSY breaking parameters are 
used between the GUT scale Mqut and a common SUSY 
threshold Msusy — ("^ti'^T-ta)^^^ (^1,2 are the stop mass 
eigenstates) , which is determined consistently with the 
SUSY spectrum. At Msusy, we impose the conditions 
for radiative electroweak symmetry breaking, calculate 
the SUSY spectrum employing the publicly available 
code SOFTSUSY JBg], and include the SUSY corrections 
to the 6-quark and r-lepton masses [lB| . The corrections 
to ttLt (almost 4%) lead (isl . IT9| to a small decrease of 
tan p. The running of the Yukawa and gauge coupling 
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constants from Msusy to Mz is continued using the SM 
rcnormahzation group equations. 

The pole mass of the top quark is fixed at its cen- 
tral value Mt = 173 GeV [33|, which corresponds to the 
running mass mt{mt) = 164.6 GeV. We adopt also the 

central value [38|] of the MS 6-quark mass mt, {mt,)^^ — 
4.19 GeV, which is evolved up to Mz using the central 
value as{Mz) = 0.1184 [sl] of the strong fine structure 
constant and then converted [s^ to the 6-quark mass in 
the DR scheme at Mz yielding mb{Mz) = 2.84 GeV. 
Finally, the tau-lepton mass is taken to be mr{Mz) = 
1.748 GeV. 

The model parameters are restricted by a number of 
phenomenological and cosmological constraints, which 
are evaluated by employing the latest version of the pub- 
licly available code micrOMEGAs jSOj. We now briefly dis- 
cuss these requirements paying special attention to those 
which are most relevant to our investigation. 

a. Cold Dark Matter Considerations. The 95% c.l. 
range for the CDM abundance, according to the results 
of WMAP [I, is 



ncnuh'^ = 0.1126 ±0.0072. 



(16) 



In the CMSSM, the LSP can be the lightest neutralino x 
and naturally arises as a CDM candidate. The require- 
ment that its relic abundance JIlsp^^ does not exceed 
the 95% c.l. upper bound derived from Eq. (1161) . i.e. 



f^LSp/i' < 0.12, 



(17) 



strongly restricts the parameter space of the model, since 
^Lsph^ generally increases with the mass of the LSP 
tolsp and so an upper bound on mLSp can be derived 
from Eq. (|17p . The lower bound on JIlsp^^ is not taken 
into account in our analysis since other production mech- 
anism s l4ll o f LSPs may be present too and/or other par- 
ticles [42|,|4J| may also contribute to the CDM. We calcu- 
late f^LSP^^ using the micrOMEGAs code, which includes 
accurately thermally averaged exact tree-level cross sec- 
tions of all t he ( co)annihilation processes [1, 13, treats 
poles [13, [3 [45] properly, and uses one-loop QCD and 
SUSY QCD corrected Higgs decay widths and 

couplings to fcrmions. 

b. The Higgs Boson Mass. According to recent in- 
dependent announcements from the ATLAS ,1] and the 
CMS l2| experimental teams at the LHC - see also 
Ref. 13 - a discovered particle, whose behavior so far has 
been consistent with the SM-like Higgs boson, has mass 
around 125— 126 GeV. More precisely the reported mass 
is (126.0 ± 0.4 ± 0.4) GeV [ij or (125.3 ± 0.4 ± 0.5) GeV 
Q. In the absence of an official combination of these 
results and allowing for a theoretical uncertainty of 
±1.5 GeV, we construct a 2 — cr range adding in quadra- 
ture the various experimental and theoretical uncertain- 
ties and taking the upper [lower] bound from the ATLAS 
[CMS] results: 



This restriction is applied to the mass ruh of the CP-even 
Higgs boson h of MSSM. The calculation of m/,. in the 
package SDFTSUSY includes the full one-loop SUSY 
corrections and some zero-momentum two-loop correc- 
tions 47]. The results are well tested [!§] against other 
spectrum calculators. 

c. B-Physics Constraints. We also consider the fol- 
lowing constraints originating from _B-meson physics: 



BR B, 




— 7> fi^fi ) of the pro- 



• The branching ratio BR {Bs 
cess Bs — > /i^A*^ 01 is to be consistent with 
the 95% c.l. bound 



< 4.2 X 10 
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(19) 



which is significantly reduced relative to the pre- 
vious experimental upper bound 32| adopted in 
Ref. |22|. This bound implies a lower bound on 
jTiLSP since BR {Bs — > ^^fi^) decreases as tolsp in- 
creases. 

• The branching ratio BR {b — > sj) of the process 
6 S7 liel, ligj is to be compatible with the 95% 



c.l. range [21|, l50|, l51| 



2.84 X 10"^ < BR(fe -> S7) < 4.2 x 10""^. 



(20) 



Note that the SM plus the and SUSY contri- 
butions [H, |4§] to BR {b sj) initially increases 
with TTiLSp and yields a lower bound on tolsp from 
the lower bound in Eq. ((20)) - for higher values of 
'T^-LSPi it starts mildly decreasing. 

• The ratio R (B„ -> tv) of the CMSSM to the SM 
branching ratio of Bu — >■ t v [s^ . [s^ is to be con- 
fined in the 95% c.l. range |5C 



0.52 < R(B„ ^ Ti^) < 2.04 



(21) 



A lower bound on mLSP can be derived from the 
lower bound in this inequality. 

d. Muon Anomalous Magnetic Moment. The dis- 
crepancy 5a^ between the measured value of the muon 
anomalous magnetic moment and its predicted value in 
the SM can be attributed to SUSY contributions aris- 
ing from chargino-sneutrino and neutralino-smuon loops. 
The relevant calculation is based on the formulas of 
Ref. [11]. The absolute value of the result decreases as 
mLSP increases and its sign is positive for ^ > 0. On the 
other hand, the calculation of is not yet stabilized 
mainly because of the ambiguities in the calculation of 
the hadronic vacuum-polarization contribution. Accord- 
ing to the evaluation of this contribution in Ref. ,27j . 
there is still a discrepancy between the findings based on 
the e^e^-annihilation data and the ones based on the r- 
decay data - however, in Ref. [30|, it is claimed that this 
discrepancy can be alleviated. Taking into account the 
more reliable calculation based on the e^e~ data 



122 < TO^/GeV < 129.2. 



(18) the recent complete tenth-order QED contribution 
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and the experimental measurements (53 | of a^, we end 
up with a 2.9 — a discrepancy 

Sa^ = (24.9 ± 8.7) x 10"^°, (22) 

resulting to the following 95% c.l. range: 

7.5 X 10"^° < 5a^ < 42.3 x 10"^°. (23) 

A lower [upper] bound on mLSP can be derived from the 
upper [lower] bound in Eq. ([23]). As it turns out, only 
the upper bound on ttilsp is relevant here. Taking into 
account the aforementioned computational instabilities 
and the fact that a discrepancy at the level of about 
3 — (T cannot firmly establish a real deviation from the SM 
value, we restrict ourselves to just mentioning at which 
level Eq. (|22]) is satisfied in the parameter space allowed 
by all the other constraints - cf. Ref. [Gj]. 

IV. RESTRICTIONS ON THE SUSY 
PARAMETERS 

Imposing the requirements above, we can delineate the 
allowed parameter space. We find that the only con- 
straints which play a role are the CDM bound in Eq. (IT71) , 
the lower bound on m/i in Eq. (fTS]) . and the bound on 
BK{Bs — )■ fj,~^^~) in Eq. p^ . In the parameter space 
allowed by these requirements, all the other restrictions 
of Sec, mil are automatically satisfied with the exception 
of the lower bound on (Ja^ in Eq. (P^ . This bound will 
not be imposed here as a strict constraint on the param- 
eters of the model for the reasons explained in Sec. IIIII 
We will only discuss at which level Eq. p2|) is satisfied in 
the parameter space allowed by the other requirements. 

In Fig. m we present the overall allowed parameter 
space in the tan/3 — A0/M1/2 plane. Each point in 
this shaded space corresponds to an allowed area in the 
Ml/2 ~ iTT-o plane (see below). The lower boundary of 
the allowed parameter space in Fig. [T] originates from the 
limit on BR {Bs — > /i+/i^) in Eq. (fTg]). except its leftmost 
part which comes from the lower bound on riih in Eq. (|18|) 
or the CDM bound in Eq. p7)) . The upper boundary 
comes from the CDM bound in Eq. dTT]) . We see that 
tan/3 ranges from about 43.8 to 52. These values are 
only a little smaller than the ones obtained for exact YU 
or the monoparametric Yukawa quasi-unification condi- 
tions discussed in Refs. [isl. [20l - [23 | . This mild reduction 
of tan/3 is, however, adequate to reduce the extracted 
BR {Bs — ?> ji^ to an acceptable level compatible with 
the CDM requirement. In the allowed area of Fig. [1] the 
parameter AqIMi/2 ranges from about —3 to 0.1. We 
also find that, in this allowed area, the Higgs mass rrih 
ranges from 122 to 127.23 GeV and the LSP mass mLSP 
from about 746.5 to 1433 GeV. So we see that, although 
m/,.'s favored by LHC can be easily accommodated, the 
lightest neutralino mass is large making its direct detec- 
tion very difficult. The range of the discrepancy 5a ^ be- 
tween the measured muon anomalous magnetic moment 



T — ' — ' — ■ — I — ■ — ■ — ' — I — ■ — ' — ■ — I — ' — ' — ■ — r 




44 46 48 50 52 

tand 



FIG. 1: The overall (shaded) allowed parameter space of the 
model in the tanP — ^o/A4"i/2 plane. 



and its SM value in the allowed parameter space of Fig. [T] 
is about (0.35 — 2.76) x 10~^" (note that 5a^ decreases as 
tan/3 or Mi/ 2 increases). Therefore, Eq. (|22p is satisfied 
only at the level of 2.55 to 2.82 — a. 

In order to get a better understanding of the structure 
of the allowed parameter space and the role played by 
the various restrictions, we will now concentrate on the 
central value of tan /3 = 48 and delineate the allowed ar- 
eas in the M1/2 — mo and ttilsp — "Tih plane for various 
values of ^0/^1/2- These allowed areas are the shaded 
areas in Figs. [2] and [3l We observe that these areas are 
very thin strips. Their lower boundary corresponds to 
— 0, where Af^ — (rrif^ — m-LSp)/'7iLSP is the rela- 
tive mass splitting between the lightest stau mass eigen- 
state f2, which is the next-to-LSP, and the LSP. The area 
below this boundary is excluded because the LSP is the 
charged f2 . The upper boundary of the areas comes from 
the CDM bound in Eq. p?)) . while the left one originates 
from the limit on BR {Bs M^A^") in Eq. HI]). The up- 
per right corner of the areas coincides with the intersec- 
tion of the lines Af^ = and fli^sph^ = 0.12. We observe 
that the allowed area, starting from being just a point at 
Ao/Mi/2 slightly bigger than —0.9, gradually expands as 
Aq/Mi/2 decreases and reaches its maximal size around 
Ao/AIi/2 ~ —1.6. For smaller Aq/ AIi/2^s, it shrinks very 
quickly and disappears just after Aq/Mi/2 = —1.62. We 
find that, for tan f3 — A8, mLSP ranges from about 983 to 
1433 GeV, while rrih from about 123.7 to 125.93 GeV. 

We will now discuss the structure of the allowed ar- 
eas in Figs, m and [31 The fact that they are narrow 
strips along the lines with Af^ = indicates that the 
main mechanism which reduces f2LSP^^ below 0.12 is 
the coannihilation of f2^s and x's. Indeed, we find that 
the dominant processes are the T2T2 coannihilations to 
bb and rf contributing to the inverse of f2LSP^^ about 
(55 — 72)% and (11 — 15)% respectively. As already no- 
ticed in Ref. @, these processes are enhanced by the 
s-channel exchange of the heavy CP-even neutral Higgs 



2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 

M,,, (TeV) 



FIG. 2: The allowed (shaded) areas in the — Tno plane for tan/3 = 48 and various ylo/Afi/2's indicated on the graph. 




FIG. 3: The allowed (shaded) areas in the mLSP — rnn plane for tan/3 = 48 and various ^o/Mi/2's indicated on the graph. 



boson 7J, with mass m//, in the presence of a resonance 
(2mLSP — TOg) - for the relevant channels, see, for ex- 
ample, Ref. [4J]. In order to pinpoint the effect of the 
-ff-pole on the T2f| coannihilations, we must track its 
position relative to the line Af^ = 0. In Fig. SI the 
dashed lines correspond to the i?-pole, i.e. to A^f = 
with Ai/ = {ruH — 2mLSp)/2mLSP for tan/3 — 48 and 
various values of Aq/Mij^ as indicated. The solid lines 

(with their dotted part included) correspond to A^^ = From Fig. gl we see that, for Aq/Mi/2 = 1 and 0, 
and the vertically [horizontally] hatched regions are ex- the i/-pole line is far from the part of the Af^ = 



eluded by the bound on BR(Ss ^^^^J■^) in Eq. 
[lower bound on nih in Eq. ([T^]. We observe that, for 
Ao/Mi/2 — 1, the lower bound on M1/2 which originates 
from the lower bound on to^ in Eq. (1181) overshadows the 
one from Eq. (jl9|) . In all other cases, however, we have 
the opposite situation. This is consistent with the fact 
that for almost fixed M1/2 and mg, the Higgs mass ruh 
increases as A0/M1/2 decreases - cf. Ref. [6; 
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FIG. 4: Relative position of the Af^ = (solid) line (with the dotted part included) and Ah = (dashed) line for tan /3 = 4i 
and various j4o/Mi/2's indicated on the graphs. Vertically [horizontally] hatched regions are excluded by the bound in Eq. (fT^ 
[lower bound in Eq. p8[) ]. The dotted areas are the overall allowed areas. 



line allowed by all the other constraints without consid- 
ering the CDM bound. Consequently, in the neighbor- 
hood of this part, the effect of the i?-pole is not strong 
enough to reduce J7lsp^^ below 0.12 via T2f| coannihi- 
lations and no overall allowed area exists. On the con- 
trary, for Ao/Mi/2 = — 1, the 7?-pole line gets near the 
otherwise allowed (i.e. allowed by all the other require- 
ments without considering the CDM bound) part of the 
= line and starts affecting the neighborhood of its 
leftmost segment, where f^LSP^^ becomes smaller than 
0.12 and, thus, an overall allowed (dotted) area appears. 
For Ao/Mi/2 = —1.3, —1.5, —1.6, the H-pole line moves 
downwards and intersects the Af^ = line with the point 
of intersection moving to the right as decreases. 
This enhances i?-pole f2f2 coannihilation in the neigh- 
borhood of a bigger and bigger segment of the other- 
wise allowed part of the Af^ = line and, thus, leads to 
f^LSP^^'s below 0.12 generating an overall allowed (dot- 
ted) area. For even smaller Ao/Mi/2's, the iJ-pole line 
keeps moving downwards and gets away from most of 
the otherwise allowed part of the Af^ — line. Also, the 
intersection of these two lines moves to higher values of 



Afi/2 and mg and the effect of the iJ-pole is weakened 
even around this intersection. So the overall allowed area 
quickly disappears as Aq/AIi^2 moves below —1.6. 

As we have seen, in the allowed parameter space of our 
model, is very close to zero. So we can restrict our- 
selves to Afj = without much loss. Note, by the way, 
that this choice ensures the maximal possible reduction 
of ^Ilsp^^ due to f2f2 coannihilations and, thus, leads 
to the maximal allowed M1/2 [or tolsp] for given tan/3, 
Aq/Mi/2, and mo [or ruh]- In Fig. [SJ we present the al- 
lowed areas in the M1/2 — ^0/^^1/2 plane for Af^ = 
and various values of tan /3 indicated on the graph. They 
are the horizontally hatched regions. Their right bound- 
aries correspond to JIlsp/i^ = 0.12, while the left ones 
saturate the bound on BR {Bs /i+/i^) in Eq. (fT9|) - cf. 
Fig. [T] The almost horizontal upper boundaries corre- 
spond to the sudden shrinking of the allowed areas which, 
as already discussed, is due to the weakening of the H- 
pole effect as A0/M1/2 drops below a certain value for 
each tan/3. The lower left boundary of the areas for 
tan /? = 44, 45, and 46 comes for the lower bound on 
in Eq. (llSp , while the somewhat curved almost horizontal 
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M,/2 (TeV) 

FIG. 5: Allowed regions in the M1/2 — j4o/Mi/2 plane for 
= and various tan /3's indicated on the graph. The dot- 
dashed lines from top to bottom correspond to = 126.5, 
126, 125, 124.5 GeV. 

part of the lower boundary of the area for tan (3 — 44 orig- 
inates from the CDM bound in Eq. (Ull). The dot-dashed 
hues from top to bottom correspond to = 126.5, 126, 
125, 124.5 GeV. We see that the rrihS which are favored 
by LHC can be readily obtained in our model for the 
higher allowed values of tan/3. 

In Table [TTl we list the input and the output param- 
eters of the present model, the masses in TeV of the 
SUSY particles - gauginos x, X2,^X3, X4> xt xf, 9, 
squarks ti, ^2, ^1, &2, wl, ur, d^, d^., and sleptons 
fi, T2,i>r, bl, e_R, i>e " and Higgses {h, H, H^, A) 
and the values of the various low energy observables 
in four characteristic cases. Note that we consider the 
squarks and sleptons of the two first generations as de- 
generate. From the values of the various observable quan- 
tities it is easy to verify that all the relevant constraints 
are met. In the low energy observables, we included 
the spin-independent (SI) and spin-dependent (SD) light- 
est neutralino-proton (x — p) scattering cross sections 
cr|p and cr|p , respectively, using central values for the 
hadronic inputs - for the details of the calculation, see 
Ref. [21[. We see that these cross sections are well be- 
low not only the present experimental upper bounds, but 
even the projected sensitivity of all planned future exper- 
iments. So the allowed parameter space of our model will 
not be accessible to the planned CDM direct detection 
experiments based on neutralino-proton scattering. We 
also notice that the sparticles turn out to be very heavy, 
which makes their discovery a very difficult task. 

In the overall allowed parameter space of our model 
in Fig. [TJ we find the following ranges for the ratios of 
the asymptotic third generation Yukawa coupling con- 
stants: ht/hr ~ 0.98 - 1.29, ht/hr ~ 0.60 - 0.65, and 
ht/hf, ~ 1.62 — 2.00. We observe that, although exact 
YU is broken, these ratios remain close to unity. They 
can generally be obtained by natural values of the real 



TABLE II: Input and output parameters, masses of the spar- 
ticles and Higgses and values of the low energy observables of 
our model in four cases (recall that 1 pb ~ 2.6 x 10~^ GeV~^). 



Input Parameters 



tan P 


48 


49 


50 


51 


-Ao/Mif2 


1.4 


1.6 


2 


2.5 


Ml/2 /TeV 


2.27 


2.411 


2.824 


2.808 


mo /TeV 


1.92 


2.295 


3.156 


3.747 


Output Parameters 


ht/h-riMcuT) 


1.117 


1.079 


1.038 


1.008 


hb/h-riMoVT) 


0.623 


0.618 


0.613 


0.607 


ht/hi{MGVT) 


1.792 


1.745 


1.693 


1.660 


fi/TeV 


2.78 


3.092 


3.823 


4.129 


Af,(%) 


1.43 


0.93 


0.1 


0.17 


AHi%) 


3.08 


1.30 


0.11 


1.76 


Masses in 


TeV of Sparticles and Higgses 




X 


1.023 


1.110 


1.309 


1.303 


xl 


1.952 


2.117 


2.489 


2.481 




2.782 


3.088 


3.815 


4.114 


xl 


2.785 


3.091 


3.817 


4.116 


xt 


1.985 


2.117 


2.489 


2.481 


xt 


2.785 


3.091 


3.817 


4.116 


9 


4.809 


5.190 


6.042 


6.040 


ii 


3.806 


4.097 


4.761 


4.781 


i% 


3.226 


3.458 


3.967 


3.902 


bi 


3.838 


4.141 


4.853 


4.947 


62 


3.763 


4.058 


4.733 


4.757 


Ul 


4.687 


5.138 


6.186 


6.483 


UR 


4.485 


4.923 


5.946 


6.257 


dh 


4.687 


5.138 


6.187 


6.483 


dn 


4.459 


4.896 


5.914 


6.227 


fi 


2.082 


2.347 


2.979 


3.293 


Tl 


1.037 


1.121 


1.310 


1.305 


Vr 


2.075 


2.342 


2.975 


3.289 


CL 


2.453 


2.819 


3.690 


4.201 


en 


2.112 


2.476 


3.339 


3.901 




2.451 


2.818 


3.689 


4.200 


h 


0.1245 


0.125 


0.126 


0.1265 


H 


2.109 


2.249 


2.621 


2.652 




2.111 


2.251 


2.623 


2.654 


A 


2.110 


2.25 


2.622 


2.652 


Low Energy 


Observables 






lO^BR (6 ^ 57) 


3.25 


3.25 


3.26 


3.26 


lO^BR {Bs fi'^fi' 


4.17 


4.15 


3.98 


4.17 


R{Bu Tf) 


0.975 


0.977 


0.982 


0.982 


10i°5a^ 


1.11 


0.89 


0.57 


0.49 




0.11 


0.11 


0.11 


0.11 


o-|J,/10-'^pb 


6.17 


4.55 


2.44 


1.75 


cr|°/10"''pb 


1.69 


1.08 


0.43 


0.28 
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FIG. 6: The complex parameters ai and 02 for various p's indicated on the graphs for the case in the second column of Table HTl 



and positive parameter p and the complex parameters ai, 
a2, which enter the Yukawa quasi- unification conditions 
in Eq. (ITSt . 

In order to see this, we take as a characteristic exam- 
ple the second out of the four cases presented in Table lU 
which yields rrih = 125 GeV favored by the LHC. In this 
case, where hi,/hr = 0.618 and ht/hr = 1.079, we solve 
Eq. psp w.r.t. the complex parameters ai, a2 for various 
values of the real and positive parameter p. Needless to 
say that one can find infinitely many solutions since we 
have only two equations and five real unknowns. Some 
of these solutions are shown Fig. [HI Note that the equa- 
tion for hb/h-r depends only on the combination pai and, 
thus, its solutions are expected to lie on a certain curve 
in the complex plane of this combination. Consequently, 
in the ai complex plane, the solutions should be dis- 
tributed on a set of similar curves corresponding to the 
various values of p. This is indeed the case as one can see 
from the left panel of Fig. IB] For each ai and p in this 
panel, we then solve the equation for ht/hr to find a2- In 
the right panel of Fig. [51 we show several such solutions. 
Observe that the equation for ht/hr depends separately 
on a2 and p and, thus, its solutions do not follow any 
specific pattern in the a2 complex plane. Note that each 
point in the ai plane generally corresponds to more than 
one points in the 02 plane. We scanned the range of p 
from 0.3 to 3 and found solutions only for the lower val- 
ues of this parameter (up to about 0.6). The solutions 
found for ai and a2 are also limited in certain natural re- 
gions of the corresponding complex planes. The picture 
is very similar to the one just described for all the pos- 
sible values of the ratios of the third generation Yukawa 
coupling constants encountered in our investigation. So 
we conclude that these ratios can be readily obtained by 
a multitude of natural choices of the parameters p, ai, 
and a2 everywhere in the overall allowed parameter space 
of the model. 



V. CONCLUSIONS 



We performed an analysis of the parameter space of 
the CMSSM with p > supplemented by a general- 
ized asymptotic Yukawa coupling quasi-unification con- 
dition, which is implied by the SUSY GUT constructed 
in Ref. jl^ and allows an experimentally acceptable b- 
quark mass. We imposed a number of cosmological and 
phenomenological constraints which originate from the 
CDM abundance in the universe, B physics {b — 57, 
Bg — > p^ p^ , and i?„ — > tv), and the mass rrih of the 
lightest neutral CP-even Higgs boson. We found that, 
in contrast to previous results based on a more restric- 
tive Yukawa quasi-unification condition, the lightest neu- 
tralino can act as a CDM candidate in a relatively wide 
range of parameters. In particular, the upper bound 
from CDM considerations on the lightest neutralino relic 
abundance, which is drastically reduced mainly by H- 
pole enhanced stau-antistau coannihilation processes, is 
compatible with the recent data on the branching ratio 
of Bs — > fJ-^ 1^^ in this range of parameters. Also, values 
of nih — (125—126) GeV, which are favored by LHC, can 
be easily accommodated. The mass of the lightest neu- 
tralino, though, comes out to be large (~ 1 TeV), which 
makes its direct detectability very difficult and the spar- 
ticle spectrum very heavy. 

The fact that, in our model, M1/2, mo, and p generally 
turn out to be of the order of a few TeV puts under 
some stress the naturalness of the radiative electroweak 
symmetry breaking. This is, though, a general problem 
of the CMSSM especially in view of the recent data on 
the branching ratio of Bg — > p^p^ and the Higgs mass 
m/i as noted in Ref. Q. Indeed, these data not only 
strongly restrict the parameter space of the CMSSM so 
as to yield very heavy sparticle masses, but also make 
the electroweak symmetry breaking less natural. 
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